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We analyze the effect of decoherence on the violation of the Clauser-Horne (CH) inequality for 
the full electron counting statistics in a mesoscopic multiterminal conductor. Our setup consists of 
an entangler that emits a flux of entangled electrons into two conductors characterized by a scat- 
tering matrix and subject to decoherence. Loss of phase memory is modeled phenomenologically by 
introducing fictitious extra leads. The outgoing electrons are detected using spin-sensitive electron 
counters. Given a certain average number of incoming entangled electrons, the CH inequality is 
evaluated as a function of the numbers of detected particles and on the various quantities character- 
izing the scattering matrix. When decoherence is turned on, we show that the amount of violation 
of the CH inequality is effectively reduced. Interestingly we find that, by adjusting the parameters 
of the system, there exists a protected region of Q values for which violation holds for arbitrary 
strong decoherence. 
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I. INTRODUCTION 

Entanglement Q is probably the most important re- 
source for the implementation of quantum computation 
and quantum communication protocols 0- Since re- 
cently, most of the work on entanglement has been car- 
ried out in optical systems using photons cavity 
QED systems 0, and ion traps Q. Solid state sys- 
tems, however, are a very attracting arena of research 
in quantum information j^, 0, @ because, in perspec- 
tive of future applications, they should allow for scal- 
ability and integration. In this light, a number of 
different realizations of entangled electrons have since 
been proposed: hybrid normal-superconducting struc- 
tures BMElllllIlQ^PiCi, superconductor 
carbon nanotubes systems [bA Il8l 19 
the Coulomb blockade regime |g, 120 
turn dots |22| . K ondo-like impurities 
bar systems [mEI I2ril27t 12 
2D electron gas |29j . 

Besides its generation, a crucial issue is that of the de- 
tection of entanglement. By means of a beam-splitter, 
entanglement can be detected in transport through an 
analysis of current noise [3(| or higher cumulants |3lj . 
Furthermore, the presence of entanglement can be re- 
vealed by analyzing the Bell inequality and quantities 
like concurrence |32(, which have been expressed in 
terms of zero-frequency charge and spin-current noise 
E3 IH El El S EH Violation of a Bell inequality 
implies that there exist quantum correlations between the 
detected particles that cannot be described by any local 
hidden variables theory. In the same spirit as it was done 
for the noise, in Ref. |35| a Clauser-Horne(CH) inequal- 
ity |32i was derived for the Full Counting Statistics 



quantum dots in 
2lJ, chaotic quan- 
, quantum Hall 
Coulomb scattering in 



(FCS) of electrons and its properties were discussed |39l |. 
In particular, it was found that the maximum violation 
of the CH inequality for electrons in the Bell state simply 
scales as the inverse of the number of injected particles. 
It was also found that the CH inequality is violated for 
a superconducting hybrid structure and, more interest- 
ingly, for a three terminal fully normal device. 

In real systems electrons are unavoidably coupled to 
the electromagnetic environment. As a result dephasing 
takes place, thereby reducing and eventually destroying 
entanglement. Understanding the consequences of de- 
phasing is an important issue. In Refs. 0, fTEl I2M liol 
the effect of dephasing was mimicked by introducing 
in the density matrix of the electronic entangled states 
a phcnomcnological parameter which suppresses its off- 
diagonal elements. By properly choosing the transmis- 
sion probability of beam-splitters or tunnel barriers, vi- 
olation of Bell inequality was found even for "strong" 
dephasing. In Refs. |2J, |4l] dephasing was introduced 
averaging over an uniform distribution of random phase 
factors accumulated in each edge channel of the quan- 
tum Hall bar. If the two edge channels are mixed by the 
tunnel barrier, no violation was reported for "strong" 
dephasing. The effect of decoherence and relaxation has 
also been analyzed using a Bloch equation formalism in 
Ref. E3 

In the present work we analyze the CH inequality for 
the FCS [35| in the presence of dephasing. We consider 
the prototype setup depicted in Fig. ^ consisting of a 
generic entangler connected to two conducting wires. En- 
tangled electrons injected in the two leads are detected 
by performing spin-selective counting along a given local 
quantization axis. The entangled electrons are subject 
to decoherence while transversing the conductors (thus 



2 



before reaching the detectors) [3(| . Various phenomeno- 
logical methods have been developed to treat dephasing 
in transport through mesoscopic conductors. In Refs. 
I43I liil which actually describes exactly nonequilibrium 
radiation acting on the system, dephasing is induced by 
a classical fluctuating potential. In Ref. |4£l dephasing is 
treated as random fluctuations of the phase of propagat- 
ing modes through the conductor. Both methods have 
been recently applied to FCS in Refs. 0, ^3 I n this 
paper decoherence is introduced as due to the presence 
of additional fictitious reservoirs along both wires. This 
method, which mimics the effect of inelastic processes, 
was introduced by Biittiker 0, |Z^| in terms of fictitious 
extra leads [$3 . The advantage of this model resides in 
the fact that inelastic, phase randomization processes are 
implemented within an elastic, time-independent scatter- 
ing problem. In the rest of the paper we shall refer to 
decoherence as to the effect produced by such fictitious 
additional leads. 

As expected, we find that decoherence suppresses the 
violation of the CH inequality, though leaving unchanged 
the range of angles for which violation occurs. In partic- 
ular, the value of the maximum violation is suppressed 
more rapidly as compared with the absence of decoher- 
ence (exponentially with the square root of the number of 
injected electrons instead of algebraically). Importantly, 
by studying the CH inequality as a function of the num- 
ber of transmitted electrons, there exist values of such 
quantity that are more protected against decoherence. 

The paper is organized as follows: In Section [H] we 
described in detail the mesoscopic system we are consid- 
ering to test the violation of the CH inequality together 
with the phcnomenological model of decoherence. Sec- 
tion [ffi] is devoted to the formulation of the CH inequal- 
ity for the FCS within the scattering approach and to 
the analysis of the no-enhancement assumption (Section 
IhTA"|) . The results are presented in Section llVl where a 
systematic analysis of the violation of the CH inequality 
against all the parameters of the device is addressed. A 
concluding summary is provided in Section Ivl 

For completeness, we include in Appendix ^ the re- 
sults relative to an asymmetric setup, whereby decoher- 
ence occurs only in one of the two wires. In Appendix 
151 and ICl we collect, respectively, the expressions of the 
expectation values and the different probability distribu- 
tions. 



II. DESCRIPTION OF THE SYSTEM 

We consider the setup illustrated in Fig. It con- 
sists of an entangler, two conducting wires and two spin- 
selective counters. The entangler, on the left-hand-side, 
is a device that produces pairs of electrons, with energy 
E < /xl, in a maximally spin entangled state (Bell state). 
On the right-hand-side of Fig. the electron counting is 
performed in leads 1 and 2 (at equal electrochemical po- 
tential /ir) for electrons with spin aligned along the local 
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FIG. 1: (Color online) Idealized setup for testing the CH 
inequality for electrons in a solid-state environment in the 
presence of decoherence. It consists of three parts: An en- 
tangler that produces pairs of spin-entangled electrons exit- 
ing from terminals 3 and 4. Two conductors that connect 
these terminals with exiting leads 1 and 2, and two analyzers. 
The conductors are described by the elastic scattering matri- 
ces Sh u , Sru, Su and Sri, and the inelastic ones Sp u and Spi. 
These last ones can simulate phenomenologically the presence 
of decoherence through the coupling via two leads (5, 6 and 
7, 8) to two additional reservoirs of chemical potentials /i and 
p! . Electron counting is performed in leads 1 and 2. Finally, 
8i and 82 are the angles at which the spin-quantization axis 
are oriented. 



spin-quantization axis at angles 9\ and 62 (spin-selective 
counters). As a convention we say that the analyzer is 
not present when the electron counting is spin-insensitive 
(electrons are counted irrespective of their spin direc- 
tion). Since we assume no back-scattering from coun- 
ters to the entangler, the particles which are not counted 
are lost and hence there is no communication between 
the two detectors. Leads 3 and 4 of the entangler are 
connected to exit leads 1 and 2 through two conductors, 
where inelastic processes are introduced through the ficti- 
tious lead model of Biittiker [^3] ■ Let us analyze in detail 
the upper wire (see Fig. which connects the emitting 
lead 3 with the exiting lead 1. The conductor consists of 
three scattering regions. The elastic scatterer connecting 
lead 3 to 3' is described by the matrix 
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[The index L (R) stands for Left (Right) elastic scat- 
terer, while P stands for Probe scatterer; u (1) refers to 
the upper (lower) wire.] Here r a (t a ) is the probabil- 
ity amplitude for an incoming particle with spin a from 
lead 3 to be reflected (transmitted into lead 3'). For a 
normal single-channel wire we set t a = t' a = ^/Tq, and 
r a = r' a = — T , where T is the transmission prob- 
ability. Inelastic scattering is introduced by plugging in 
an additional reservoir of chemical potential fi with an 
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energy- and spin-independent scattering matrix 
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represented by a triangle in Fig. ^ For the sake of clar- 
ity, we have denoted by a check ("), a caret (~) and an 
overbar (~), respectively, (2x2), (4x4) and (8x8) matri- 
ces. In Eq. 1 and 6 are, respectively, unit and zero 
(2x2)-matrices in the spin space, and (1 — a) is the prob- 
ability for transmitting a particle between leads 3' and 
1'. The coupling parameter ranges from a = 0, when no 
particles are transmitted into leads 5 and 6 from leads 3' 
and 1', to a = 1, when no particles are transferred be- 
tween leads 3' and 1'. A third elastic scatterer, described 
by a matrix Sr. u , connects lead 1' to lead 1. The con- 
ductor is therefore described by the matrix S13 defined 
as S 13 = Sp n ® Sp u <g> Sl u , where the notation g) stands 
for the scattering matrix composition (elimination of in- 
ternal current amplitudes) |53(. For simplicity, we shall 
assume that Sr u = Slu- 

Due to the presence of the additional reservoir, parti- 
cles propagating through lead 3 are transmitted partially 
to lead 1 and partially to leads 5 and 6 (see Fig. ^) . The 
additional reservoir, however, can transfer itself particles 
to lead 1 and 3. As a result, only a fraction of the parti- 
cles arriving in lead 1 comes from coherently transmitted 
ones sent in from lead 3, with probability 



Tl3 = 



T 2 (l-a) 



[l + (l-T )(l-a)P 



(3) 



Another fraction, the incoherent contribution, comes 
from the additional reservoir through leads 5 and 6, with 
probability 



Tib + T 



T a 
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l + (l-T )(l-a) 



(4) 



The presence of the extra reservoir mimics the fact that 
the current flowing through the conductor is partially 
composed of particles (the incoherent fraction) which 
have lost phase memory while traversing it. For a = all 
particles are coherently transmitted and T 15 + T 16 = 0, 
while for a = 1 all particles are transferred incoher- 
ently and T13 = 0. For a = 0, the overall trans- 
mission probability through the conductors is given by 
T = T 2 /(2 - T ) 2 . In the rest of the paper we will refer 
to a as to the decoherence rate. 

The chemical potential /1 of the additional reservoir is 
set in such a way that no net current flows in or out of 
the reservoir (I5 + Iq = 0). This constraint is enforced 
only on average. An instantaneous current in or out the 
additional reservoir is then allowed and 
a non-fluctuating chemical potential fi is assumed (for 
this reason the additional terminal does not behave as a 
voltage probe). 



A similar description applies to the lower wire connect- 
ing lead 4 with lead 2, so that the scattering matrix of 
the conductor is defined as §24 = Sri (g> Spi <g> Su , where, 
for simplicity, we set Sri = Su- If the angles 6\ and 
#2 of the analyzers are parallel to each other and in the 
absence of spin mixing processes, the total matrix of the 
system can be written as 



S 



Sl 3 

6 s 24 



(5) 



The general scattering matrix relative to non-collincar 
angles Se 1 .e 2 is obtained from S by rotating the spin 
quantization axis independently in the two conductors 
(note that this is possible because the two wires are de- 
coupled) j5|| : Sg 1: g 2 = USW , where U is the rotation 
matrix given by 



U 



1 1 6 \ 

b u 6l 
6010 

V U 02 J 



(6) 



and 



Un 



cos ■ 



sin ■ 



(7) 



For simplicity, we further assume that the two conduc- 
tors are equal and that they are subjected to the same 
degree of decoherence, so that S13 = S24. For this reason 
the chemical potentials of the additional reservoirs are 
identical. It is interesting to notice that decoherence pro- 
cesses in the two wires are, in some sense, " uncorrclatcd" , 
meaning that we have imposed that the currents flowing 
through the fictitious leads vanish separately in the two 
reservoirs. (Correlations can be introduced, for example, 
by imposing ^5 + I§ + 17 + 1$ = 0.) In the symmetrical 
setup we are considering here, fj, = /xr, + (ml — A*r)/2. In 
the rest of the paper we consider the case in which all 
the reservoirs are at zero temperature. 

The incoming state of the system | ip) depends on 
whether the energy of electrons falls within the range 
Mr < E < or /x < E < 



IV>) 



M < E < ^ L 
Mr < E < fi 



(8) 



where 



4 1 (E)al l (E)±al l (E)al,(E)] 0) , (9) 



and 



ife) = ■^=[al 1 (E)al l (E)±al l (E)al 1 (E) 



n=5,6,7,8 



(10) 
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In Eqs. © and l|10|) a\ a (E) is the creation operator for 
a propagating electron in lead i with spin a at energy E. 
The upper sign refers to the case in which the incoming 
state is a spin triplet and the lower sign to the spin sin- 
glet. Electrons with energy between fi and fi^ are exiting 
leads 3 and 4 of the entangler in a superposition of spin 
| and I states. For energies between /^r and fi electrons 
are also injected from the additional leads (with indexes 
5, 6, 7 and 8) in a factorized state. Note that this occurs 
only in the presence of decoherence, i.e. for q^O. 

By setting /xr = and fi^ = eV, the total current 
flowing in lead 1, calculated using the Landaucr-Buttikcr 
formalism in the linear response regime, is given by 



h = e 2 V/h(T 13 + T 15 + T] 



IB I 



(11) 



Although the coherent part of the current decreases with 
a, the total current increases with it (except for T = 
1, where it remains constant and equal to e 2 V/h). We 
would like to mention that this is a special feature of the 
model of decoherence we are using, not to be expected in 
general. 



III. CH INEQUALITY FOR THE FULL 
COUNTING STATISTICS 



The quantity employed in the formulation of the CH 
inequality, as derived in Ref. l35l is the joint probability 
P(Ql, Q 2 ) for transferring a number of Q\ and Q 2 elec- 
tronic charges into leads 1 and 2 over an observation time 
t. The CH inequality is based on the hypothesis that the 
outcome of a measurement could be accounted for by a lo- 
cal hidden variable theory. The test of the CH inequality 
proceeds as follows. The entangler is switched on during 
an observation time t (where the minimum t is the inverse 
of the measuring device bandwidth) in which it emits an 
average number M of pairs of entangled electrons. After 
traversing the conductors (and being affected by inelastic 
scattering) the electrons are counted in both terminals 1 
and 2. The experiment is then repeated to get single ter- 
minal and joint terminal probability distributions that 
Qi particles arrive into analyzer 1 and Q 2 particles ar- 
rive into analyzer 2 (along a local spin-quantization axis 
or independently of it) with Qi + Q 2 < 2M. 

The CH inequality for the FCS reads 

s CH = P 9l ' 92 (Qi,Q2) - P 9l ' 9 ' 2 (Qi,Q2) + P e ' 1 ' e '(QuQ2) 

+P 9 '^(Q 1 ,Q 2 ) - P e '^-(Q 1 ,Q 2 ) 
-P-'" 2 (Qi,Q 2 )<0. (12) 

The possible violation, or the extent of it, also depends 
on Qi and Q 2 . P Bl,e2 {Qx,Q 2 ) is the joint probability 
in the presence of two analyzers, where Q± electrons are 
counted in lead 1 along B\ direction and Q 2 are counted 
in lead 2 along 9 2 . P 01 <~ (Q 1; Q 2 ) is the corresponding- 
joint probability when one of the two analyzers has been 
removed. The same notation will be used for single termi- 
nal probability distributions: P 0i (Qi) in the presence of 



an analyzer and P(Qi) if no analyzer is present. Eq. i|12|> 
holds for all values of Qi and Q 2 which satisfy the no- 
enhancement assumption: 



P 9 ^Qi)<P{Q l ). 



(13) 



The joint probability distribution for transferring Qi a 
electrons with spin a in lead 1, Q 2a electrons with spin 
a in lead 2 and so on is given by 



P(Qi h Q n ,Q 2h ...) 

1 r+ 7T 



(27T) 



2n 



2i 



x(At,Aj.) 



(14) 



where x(A|, AjJ is its characteristic function that can be 
expressed within the scattering approach. 

For long measurement times t, the total characteristic 
function \ is the product of contributions from different 
energies, so that 



(15) 



The energy-resolved characteristic function for the trans- 
fer of particles at a given ener gy E in a structure attached 
to n leads can be written as 153, |59, 1601 



xs(At, Aj.; 



n 



x ne 



(16) 



where the brackets (...) stand for the quantum statisti- 
cal average over the thermal distributions in the leads. 
Assuming a single channel per lead, ^0(1) lii ^ ne num b er 
operator for outgoing (incoming) particles with spin a in 

lead j and Af , A^ are vectors of n real numbers, one for 
each open channel. In terms of outgoing (incoming) cre- 
ation operator (aj CT ), which are linked by the total 
scattering matrix of the system S, the number operators 
can be expressed as 



(17) 



At zero temperature, the statistical average over the 
Fermi distribution function in Eq. i|16|) simplifies to the 
expectation value calculated over the state | ip) defined 
in Eq. JSJ. The interval of integration in Eq. I|15|) can 
be separated in two energy ranges, namely E < fi and 
fi < E < eV. Since, in the limit of a small voltage bias V, 
Xe is energy-independent, Eq. 1|15|) can be approximated 
to 



x(^T. A l) 



M„ r _ _ -1 M-Mu. 

X g(A T ,A x ) Xo B (At,A x ) , (18) 



where A/ M = fit/h and M — eVt/h. 
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According to Eq. (|14f> . both single terminal and joint 
probability distributions require the computation of mul- 
tidimensional integrals, which can only be performed nu- 
merically. In Appendix |C] it is shown that the various 
probability distributions needed to evaluate the CH in- 
equality can be expressed in a differential form, more 
suitable for numerical evaluation. All the expectation 
values needed for the calculations are collected in Ap- 
pendix [B] Since the two wires are decoupled and there 
are no spin-flip processes, the joint probabilities with a 
single analyzer are factorized: 

P 9l >-(Qi,Q2) = P 6l (Qi)P(Q2), 
p-- 92 (Qi,Q2) = p(Qi)p 02 (Q2). (19) 

Rotational invariancc makes P 8l -~(Qi, Q 2 ) and 
P~' 02 {Qi,Q 2 ) independent of the angle of the analyzers, 
while P 01 ' 02 (Qi,Q 2 ) depends on the angles only through 
the combination ei ^ 8 ' 2 (upper sign for triplet and lower 

e 1 ±e 2 

sign for singlet), so that we can define P 1 2 (Qi, Q2) = 
P 01 ' 02 {QuQ2) and P 1 ^(Q 1 ,Q 2 ) = P 6 ^-'(Q 1 ,Q 2 ). Asa 
result, the CH inequality depends only on three angles 
a = 0i±0 2 , Ob = 02 ±#1 and C = 6\±6' 2 (0 d = 0i±0 2 is 
a linear combination of the other three: 9d = a + 0b + c ). 
Since P 9l ' 02 (Qi, Q 2 ) is an even function of in or- 

der to find maximal violations we can restrict the evalu- 
ation of the CH inequality to the following set of angles: 
6a = 05 = C = 0d/3 = 29. (This is found by im- 
posing that positive contributions to 5ch are maximum 
while negative contributions are minimum.) The quan- 
tity 5ch , characterizing the CH inequality, will therefore 
depend on a single angle 0, on the decoherence strength 
a and on the value of the transmitted charge Q\ and Q 2 . 
As a result, the CH inequality takes the simplified form 

Sen = 3P 1 e 2 (Qi,Q 2 )-P 1 3 |(Q 1 ,Q 2 )-P 1 ,_(Qi,Q 2 ) 
-P_, 2 (Qi,Q 2 ) <0. (20) 

Without loss of generality we can choose a = a' =1 . The 
other cases can be recovered by rotating the polarizers 
an angle tt. 

A. No-enhancement assumption 

As mentioned above, the CH inequality can be de- 
rived under the no-enhancement assumption, Eq. (|f 3fl . 
Such a condition is trivially true when a single particle 
is transmitted, Q = 1: The presence of an analyzer can 
only decrease the counting probability |37| . However, 
when many particles are transmitted, Q > I, the no- 
cnhanccmcnt assumption is a relationship between dis- 
tribution probabilities that is not, in general, satisfied 
for all values of Q. 

We remind that in the absence of decoherence [35|, 
for given M and Q, the no-enhancement assumption in 
one of the two leads is satisfied only within a range of 
values of T below certain threshold T max (M, Q). In the 




FIG. 2: (Color online) a) Maximum value of the transmis- 
sion, T m ax, allowed by the no-enhancement assumption as a 
function of decoherence rate a for M = 20 emitted pairs and 
for the different values of Q. b) Minimum allowed number of 
transmitted particles Q m i n for a fixed wire transmission as a 
function of the decoherence rate. 

case of different numbers of transmitted particles in lead 
1 and 2 (Qi ^ Q 2 ), the maximum allowed transmission 
probability must be taken to be the minimum between 
T ma .x(M, Qi) and T max (M, Q 2 ), according to our assump- 
tion of identical wires. 

The no-enhancement assumption is affected by deco- 
herence as a consequence of the fact that single termi- 
nal probabilities, with or without analyzer, depend on a. 
More precisely, the no-enhancement assumption in one of 
the two leads is satisfied for transmissions up to a thresh- 
old value which is now a function of a: r max (a, AI, Q). 
Unlike the ideal case, for a 7^ it is not possible to find 
an analytical expression for T max . In Fig.JJJi P m ax is plot- 
ted as a function of a for M = 20 and all values of Q 
from 1 to 20. One can see that r ma x monotonically de- 
creases with a for values of Q < M/2 and monotonically 
increases for large values of Q. For intermediate values of 
Q, Pmax decreases up to values of a close to one and then 
rapidly increases reaching one when a = 1. This behavior 
is specific of the fictitious lead model and reflects the fact 
that both the average total current [Eq. (|ll|l ] , related to 
P(Q), and the average spin-polarized current, related to 
P e (Q), are increasing functions of a. Indeed, as a conse- 
quence of a finite a, the two distributions shift to larger 
values of Q, as it would happen for an enhanced effective 
transmission probability. Its maximum allowed value by 
the no enhancement assumption is therefore reached for a 
smaller T. As a consequence T max must decrease with a. 
This argument is not valid when T max ~ 1 at a = 0, since 
the average currents do not change appreciably with a 
and only the peculiar shape of the distributions matters. 
We define 

Pmax(a, M, Qi, Q 2 ) 

= Min[T max (a, M, Qi), T max (a, M, Q 2 )]. (21) 

Alternatively, given a wire with a fixed transmission T, 
the no-enhancement assumption is verified for values of 
Q bigger than or equal to a certain value Q m i n (a, M, T). 
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FIG. 3: (Color online) The quantity Sch is plotted as a func- 
tion of O for M — 20, Qi = Qi = 8 and different decoher- 
ence rates at the corresponding maximum allowed transmis- 
sions. In particular, for a = (solid line) T max = 0.54, for 
a = 0.3 (dashed line) T max = 0.43 and for a — 0.5 (dot- 
ted line) T max = 0.35. The amount of violation of the CH 
inequality decreases with a, whereas the range of angles for 
which violation occurs does not change. We call Obost the 
angle corresponding to the maximum violation. 



For a 7^ 0, the behavior of Q m [ n is shown in Fig. |2b 
for M = 20 and for different transmissions. We observe 
that it increases (step-wise, since only integer values of 
the number of particles are permitted) as a function of 
the decoherence rate for almost every transmission T, 
except for those close to unity, for which it decreases. The 
behavior for small values of T can still be understood in 
terms of the average current increase with a. For T = 1, 
being Q m i n = M at a = 0, decoherence can only cause a 
decrease. 



IV. RESULTS 

In the present section we shall discuss how the CH 
inequality of Eq. (|20J) is affected by the presence of de- 
coherence. There are some general characteristics of the 
behavior of Sc xt th at were already found in the absence 
of decoherence [35| that hold also for finite a j6j| . The 
most relevant are the following: 

• <Sch is always symmetric function of around 
9 = tt/2; 

• for given M, Qx and Qi the maximum violation 
always occurs for T = T max (a, M, Qx, Q2)- 

In Fig. |21 Sea is plotted as a function of the angle 
for M = 20 and Q\ = Q2 — 8. The three curves refer, 
respectively, to a = (solid line), a = 0.3 (dashed line) 
and a = 0.5 (dotted line), each one calculated for the 
corresponding T = T max reported in the label box. The 
plot shows that the CH inequality is violated within a cer- 
tain window of values of 0. The violation is suppressed 
with increasing decoherence rate, but occurs for the same 
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FIG. 4: Density plots of the maximum value of Sch, eval- 
uated at T = r max (a, Qi, Q2) and = 6best(a, Qx, Q2), in 
the (Qi,Q2) plane for M = 20 relative to four different val- 
ues of decoherence (a = 0, 0.2, 0.4, 0.6). Sch ^ 0.012 is 
the maximum violation for M — 20 found in the absence of 
decoherence. 



range of angles. This is due to the following proper- 
tics of the joint probabilities, which hold at T = T max 



for all values of a: i) Pxf^iQuQ^) = Pi^^iQi, Q2) = 
Pi ~{Qi,Qi), as a consequence 5ch(0 = 7r/4,37r/4) = 
0;'and ii) P 1 e 2 (gi,Q 2 ) > iff (Qi, Q2), Pi,-(Qi,Q 2 ), 
P-AQx,Q2) for tt/4 < < 3tt/4. We checked that 
by reducing T from T max , but keeping a constant, both 
the window of angles where violation is present and its 
amount are decreased. Note that between = and 
7r/2, there is always an angle for which Sen is maximum, 
we shall denote it by 0b es t(o;, M, Q±, Q2)- For given a, 
M, Qi and Q2, the maximum violation occurs at T max 
and ©best- 

We now analyze the maximum violation of the CH in- 
equality for a given M with T = T max and = ©best 
as a function of Qx,Q2 and a. In Fig. ^we show four 
density plots of Sen in the (Qx,Q2) plane for different 
values of decoherence rate and M = 20. In the gray scale 
white corresponds to Sen = and black to its maximum 
value taken in the absence of decoherence. When a = 0, 
the CH inequality is strongly violated for diagonal terms 
of the distribution (where Q\ = Q2)- However, some 
weaker violations are also possible for Qx 7^ Q2, though 
they tend to disappear with increasing a. By increasing 
the values of a the plots show that the maximum viola- 
tion of the CH inequality decreases rapidly: for a = 0.6 
we get only 16% of the largest value reached at a = 0. 
The behavior of the CH inequality is symmetrical with 
respect to the exchange of Qx with Q2 for any rate of 
decoherence. In Fig. [S] we report the section of the plots 



FIG. 5: (Color online) Maximum value of the quantity Sen 
(T — T max and = Gbest) as a function of Q for M = 20 
and different values of decoherence: a = 0, 0.3, 0.5, 0.8. The 
largest violation always occurs in the absence of decoherence 
and, for a given Q, violation is reduced monotonically with 
a. The position where the maximum occurs, Qbcst, slightly 
increases with a. At Q = Qbcst, Tmax = 0.54 for a = 0, 
Tmax = 0.59 for a = 0.3, T max = 0.52 for a = 0.5 and T max = 
0.39 for a = 0.8. 



FIG. 6: (Color online) Sch, normalized to its value in the 
absence of decoherence and calculated at T = T ma x, O = 
Gbcst and Q = Qbcst, is plotted as a function of decoherence 
rate a and number of injected entangled pairs M. Longer 
measuring times (i.e. larger values of M) make decoherence 
more effective, that is, make the detection of entanglement 
more difficult. 



in Fig. 21 along the diagonal of the (Qi, Q2)-plane. The 
four curves are relative to a = 0, 0.3, 0.5 and 0.8 an 
M = 20. Several observations are in order. If we denote 
with Qbcst the position of the maximum of a curve, for 
all values of decoherence rate Qbest ~ M/2, more pre- 
cisely, Qbest = 8 for a = and Qbcst = 10 for all other 
curves. This slight increase of Qbcst with a is due to the 
fact that an increase in decoherence is accompanied by a 
slight enhancement of the average current [Eq. (|llfl ] flow- 
ing through the wires (as mentioned at the end of Section 
HI)l . This is, however, a specific feature of the model of 
decoherence we are considering. Note furthermore that, 
as decoherence gets stronger, the range of values of Q for 
which violation takes place shrinks. 

We now discuss the violation of the CH in- 
equality as a function of a and M at T = T max , 
© = ©best and Q = Qbcst- In Fig. the ratio s = 

SCH (a, M, T max , Obest, Qbcst )/<5>Ch(0, M, T max , Obest, Qbest) 

(i. e. the quantity Sqh normalized to its value in the ab- 
sence of decoherence) is reported in a three-dimensional 
plot as a function of a and the number of emitted 
pairs M. The most interesting feature is that such a 
ratio decays more rapidly with a as M is increased. 
This means that decoherence is more disruptive, as 
far as detection of entanglement is concerned, for long 
measuring times (i.e. large M). As an example, for 
AI = 1000 the extent of the violation is reduced by 80% 
at a = 0.1. More precisely, for values of M larger than 
30, we find that the normalized Sen follows the law: 



sinh[ J fC(l - a) bVM ] 
sinh(if) 



(22) 



with K = 7.26 and b = 0.076. 

Another interesting aspect is related to Qbcst which, 
as mentioned above, only slightly increases with a for 




FIG. 7: (Color online) Qbcst is plotted as a function of M for 
different values of decoherence rate (a = 0, 0.2, 0.5 and 0.8). 
In the inset, curves are shown over an extended range, up to 
M = 1000. For a given a, with increasing M the value of 
Qbcst increases very slowly remaining of the order of 10. 



all values of M. As M is increased, the value of Qbcst , 
for a given a, does not increase proportionally to M, 
but very much slowly and surprisingly remains of the 
order of 10 for M = 1000 (see Fig. 0. For a = this 
can be understood as follows. On the one hand, one 
expects Qbcst, corresponding to the largest iSch, to be 
about the position of the maximum of joint probability 
distributions, which can be assumed to be equal to the 
product MT. On the other hand, T max is a decreasing 
function of M, in fact it decays as 1/M [35j. The product 
MT max is therefore expected to be a constant. Indeed, 
it is possible to show, in the large M expansion, that 

MTrr 



'best 



for a = and Q be st ~ Af y/T max for 



a 7^ 0, while T n 



1/M for a = and T u 



1/M 2 
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FIG. 8: (Color online) Maximum value (a max ) of the deco- 
herence parameter for which there is still violation of the CH 
inequality as a function of Q for M = 20. The line with ♦ 
represents a max (Q) with Sth = 0: from Q — 1 to Q — 14 vi- 
olation is found for any decoherence rate. The line with ir is 
instead computed using a threshold Sth which corresponds to 
0.25% of the maximum violation value for a = 0, and the line 
with ■ using Sth corresponding to 1%. The latter threshold 
is used in the inset where a ma x(Q) is plotted for M — 10, 20, 
30 and 40. Interestingly we found that there is a value Q that 
is more robust against decoherence. In particular, Q = 6 for 
M = 10, Q = 11 for M = 20, Q = 13 for M = 30 and Q = 14 
for M — 40. With increasing M, a max diminishes slowly. 

for a / 0. As a result, Qbost is roughly constant as a 
function of M and a [62| . 

The final point we address is the maximum deco- 
herence rate, that we denote by ev max , for which there 
is still violation of the CH inequality function of 

Q = Qi = Q2- In Fig. [S] we plot a max as a function of 
Q for M = 20 at T = T max and 6 = e bos t- The line 
with the symbol ♦ shows that violation of the CH in- 
equality is found for any rate of decoherence for Q = 1 
to Q = 14 and thereafter a max decreases sharply. Never- 
theless, the extent of violation for a close to 1 is almost 
negligible for most of 1 < Q < 14. One can therefore 
introduce a small positive threshold 5th which defines 
the violation as: 5ch < <5th- The line with -A- refers to 
a threshold of 0.25% of the maximum value of Sen at 
a = (5 th = 3 x 10~ 5 for M = 20), and the line with 
■ to a 1% (5 th = 1.2 x 10~ 4 for M = 20). The latter 
percentage is used for the thresholds of the plots in the 
inset of Fig. EI where ev max is plotted for M = 10, 20, 30 
and 40. It is shown that there are values of Q that are 
more resistent to decoherence, i. e. for which violation 
survives for larger decoherence rates. In the caption of 
Fig. [21 the most protected value against decoherence is 
denoted with Q = <2best(a m ax)- 

V. CONCLUSIONS 

In this paper we have studied the effect of decoher- 
ence on the violation of the CH inequality formulated 
in terms of the FCS J35|. The system under investiga- 



tion (Fig.m consists of an idealized entangler connected, 
through a pair of identical mesoscopic wires, to spin- 
selective counters. We have assumed that decoherence, 
which occurs equally but independently in the two con- 
ductors, is produced by the presence of additional ficti- 
tious reservoirs according to the phenomcnological model 
of Biittiker 01^3- Decoherence is parameterized by the 
rate a. 

As expected, decoherence gives rise to suppression of 
the violation of the CH inequality. The extent of such a 
suppression has been analyzed as a function of the param- 
eters which characterize the system, namely the trans- 
mission T of the wires, the angle between analyzers 0, 
the number of injected entangled pairs M and the num- 
ber of transmitted particles Q\ and Q2 in the counters. 
First we have discussed the no- enhancement assumption, 
a condition that needs to be satisfied in both leads 1 and 
2 in order for the CH inequality to hold. We have found 
that such condition, in a given lead, is verified for all 
transmission T up to some maximum value T max which 
depends on Q, M and, of course, a. In particular, T max 
decreases with the decoherence rate up to some value 
of Q and thereafter increases. The main results can be 
summarized as follows: 

• Maximal violation, even in the presence of deco- 
herence, occurs at the largest allowed transmission 
T = Tmax and for Qi = Q2 (it disappears very 
rapidly when Qi ^ Q2)- 

• As long as T — T max , the angle range of the analyz- 
ers for which violation takes place does not depend 
on the decoherence rate, though the extent of vio- 
lation decreases with a. 

• In the absence of decoherence, the maximum vio- 
lation of the CH inequality was proved to decay as 
1/M j3f|. Here we have found that, for finite a, the 
parameter 5ch decreases exponentially with vM, 
more precisely as [J(a)\*** /M, i.e. decays both 
with increasing M and a. 

• The value of Q for which maximum violation occurs 
is virtually independent of M, which means that 
the largest violations appear for relatively small 
numbers of transmitted particles, even at large ob- 
servation times. 

• Interestingly, we have found that the largest de- 
coherence rate for which the CH inequality is vi- 
olated (within a given small tolerance) presents a 
maximum as a function of Q. This means that 
there exist numbers of transmitted charges which 
are more protected against decoherence, i.e. the in- 
fluence of the environment is less disruptive as far 
as the violation of CH inequality is concerned. 

Although, in this paper, dephasing is assumed to be 
produced by the presence of additional reservoirs, other 



9 



reservoir \i 




o 



Nft 



Q, 1 







































































a=0.4 










































































































































































































































































































































































































1 









































































































































































































































































Qv 







































































a=0.6 



































































































































































































































































































































































































































































































































































































































































































W CH 
6.6 10 3 



FIG. 9: (Color online) a) Idealized setup with a single ad- 
ditional reservoir in the upper branch of the system. Scat- 
tering matrices are chosen so that, in the absence of de- 
coherence, the two conductors have equal transmission, b) 
Density plot of the maximum value of Sch in the (Qi,Q2)- 
plane for M = 20 and for a = 0.4 (left) and a = 0.6 
(right). As decoherence increases, the maximum violation 
is not achieved on the diagonal, but is shifted towards the 
right-bottom part of the plane (Qi > Q2)- This occurs 
because only the current flowing through the conductor af- 
fected by decoherence is modified. Furthermore, the sup- 
pression of the violation by a is less pronounced with re- 
spect to the case with two additional reservoirs. For ex- 
ample, for a = we find Max[ScH] = 0.012 achieved at 
(Qi = 8, Qi = 8). For a — 0.2 and one additional reservoir 
we have Max[ScH] = 0.0089 reached at (8, 7), whereas for two 
additional reservoirs we get Max[ScH] = 0.0074 at (9,9). For 
a — 0.4 and one additional reservoir, Max[ScH] = 0.0066 at 
(9, 7), and with two additional reservoirs, Max[ScH] = 0.0042 
at (10, 10). Finally, for a = 0.6 and one additional reser- 
voir, Max[iScH] = 0.0046 at (10,7), and with two additional 
reservoirs, Max[S C H] = 0.0020 at (10, 10). 



different sources of decoherence are possible in meso- 
scopic systems. We believe that this model captures the 
main effects of decoherence, as far as violations of the CH 
inequality in a mesoscopic system is concerned, and that 
the results found in this work may be useful to design 
the best experimental conditions. 

Since real systems cannot be perfectly shielded from 
the environment, the issues analyzed in this work seem 
adequate not only from a fundamental point of view, but 
also in what it might contribute to the understanding of 
the properties of lossy quantum channels. For the future 
it would be interesting to apply our method to realistic 
systems, like normal or superconducting beam splitters. 
Of interest would also be the combined effects of the pres- 
ence of spin-flipping processes and decoherence. 




FIG. 10: (Color online) a) Idealized setup with no decoher- 
ence but differently transmitting upper and lower conductors, 
b) The density plot of the maximum value of the quantity Sch 
is shown in the (Qi, Q2)-plane for M = 20. 
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APPENDIX A: ASYMMETRIC SETUP: ONE 
ADDITIONAL RESERVOIR 

It is interesting to consider the case in which decoher- 
ence affects the two wires differently. In this appendix we 
study the case when decoherence affects only one of the 
two conductors, i.e. in the presence of a single additional 
probe, for example, in the upper branch, as depicted in 
Fig. [jju Being To the transmission of the elastic scatter- 
crs in the upper conductor, we choose the transmission 
of the lower conductor to be equal to T = Tq /(2 — To) 2 , 
in order for the two conductors to have the same con- 
ductance in the absence of decoherence. Fig. |5Jd shows 
the density plot of the maximum value of Sen (with 
T = T max and 6 = 6b 0S t) as a function of Qi an( l Q2, 
for M = 20, a = 0.4 (left) and a = 0.6 (right). In the 
presence of decoherence the maximum violation is not 
achieved on the diagonal (Qi = Q2), i-e. the behavior 
of Sen is not symmetrical anymore with respect to Qi 
and Q2- This is due to the fact that, as we have seen 
above, the overall current increases with a so that it is 
more likely to transmit a larger number of particles in the 
conductor subjected to decoherence. Another difference 
with respect to the case with two additional reservoirs 
is that the suppression of the violation by a is less pro- 
nounced. 

The asymmetry found in the behavior of the density 
plots of Fig. EJ) must not be confused with the asymme- 
try we would obtain in a setup without decoherence but 
with conductors of different resistance (system sketched 
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in Fig. \Wh). In this case, by applying separately the no 
enhancement assumption to the two conductors, large vi- 
olations occur in a vast region of the (Qi, Q2)-pla,ne, as 
shown in Fig. 110b . Interestingly, we note that one would 



get large violations of the CH inequality for Qi ^ Q^. 
However, this asymmetry would not come from the fact 
that entanglement is weakened by dephasing in only one 
wire. 



APPENDIX B: EXPECTATION VALUES 

The most general expression for the characteristic function when spin-cr electrons are counted in lead 1 and spin-a' 
electrons are counted in lead 2 is 

X E (X la ,\ 2a ,) = 1 + (e~ iX ^ - 1) (N^) + (e- 4A -' - l) + (e^ - l) (e^-' - l) (N^N^') (Bl) 

for each relevant energy range: < E < fx and /i < E < eV. When both spin species are counted in one of the 
terminals, the characteristic function reads 



Xe(Xi, A 2 ,0 = 1 + (e-' tAl - 1) (« T + A#)) + (e-* A -' - 1) (N*° ) + (e"* Al - 1) ( 



+ (e-^ - l) 2 (N^N^ 1 ) + (e-' iAl - 1)' ( 



l)(«+OiV^} 
(B2) 



for counting both spins in terminal 1, where we have set Aij = Xi[ = X±. 

Using Eqs. lO, CHI) and {Bljl . at zero temperature, one can calculate the single terminal probability distribution: 



p 6i (Qi 



2n 



dX 



la 



i) m a )s 



1) (N^)b 



M-M u 



a i\laQl 



(B3) 



where {Nq')s,b = (^s,b \ Nq 7 \ ips,B)- After integration over Alcr we get 



/ « , \ Qla 



M u 



Wl\n[M ^ ,Q i a ] 

V 

n=Max[0,Qi r - M+M^] 



M - M„ 
Qia - n 



(Nb°)s(l-(N??)B) 



(B4) 



If one chooses /il = eV and /jr = 0, one obtains \x = ^f- and M M = 4^. Therefore, in order for M M to be integer, M 
must be an even number. 

Using Eqs. (|14[) . (|18|l and (|B2() . the single terminal probability distribution when both spin species are counted in 
the terminal is 



P(Ql) = h f V + (e " Al _ l) { ^ o)s + (e " Al " 1)2 &oXo}s] M " 



x 1 + (e-' lAl - 1) (N&b + (e~' lAl - l) 2 (NgN# )b " e lXlQl d\ 1 



M-M u 



(B5) 



where {N& s ,b = + N t 



o))s,B- 



In the following subsections we collect all the expectation values needed to work out the expressions for the proba- 
bility distributions of Eq. (|2U|l . 
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1. Setup with two additional reservoirs 

Let us consider the setup depicted in Fig. ^ where two additional reservoirs are present. The expectation values 
for energies < E < /i as a function of transmission T, decoherence parameter a and analyzers' angle arc 

/KrV\ /at^\ 2VT[Vr + a(2 - VT) - q 2 (l - VT)] 

{N ° )s = {N ° )s = MWf ' ( } 

{N ° )S ~ {N ° )S [2-*(l-VTW ' (B7) 



[2 - a(l - VT)] 3 



>1T - 2T 4T{ [a(2 - a) + VT(1 - a + a 2 )] 2 - T(l - a) 2 cos 6} 

v v o iV o )s = ^ t: /^rr; ' ( By ) 



[2 - a(l - VT)] 4 

12 



<« = < = 8r[ ^ + ;%fj-;r^ )]2 . (bio 



(Bll) 



(N£n^N£) s = (N^N^N* 1 ) 3 = 8VT 3 "a[2Vr + 2a(l-VT)-a 2 (l- VT)] 

^ [VT + a(2 - VT) - a 2 (l - VT)] 
[2-a(l-VT)p 
The expectation values for energies [i < E < eV are 

= ^'''' p-XV (B12) 

<~°>= - ( *°>» - F^cT^W ' (B13) 

(»<!,' ffo l )B = (W?JV^) B = 0, (B14) 
V ° 07 [2-a(l-VT)] 4 

(*?JV^JVj r ) B = (Wi'wJ'ivg'jB = 0. (B17) 



2. Setup with one additional reservoir 

Let us now consider the asymmetrical setup of Fig. [51 where there is only one additional reservoir. For energies 
< E < fi we have that (Nq)s, (Nq)s and (Nq Nq}s are equal to the case with two fictitious reservoirs. The other 
expectation values are 

(N$)s = §, (B18) 

(TV 2 ) s = T, (B19) 

(N^N 2 l ) s = 0, (B20) 

= V^ Q (2- a ) + T 2 [l- Q + . 2 -(l- a) cose] 
V ° 07 [2 - a(l - VT)] 2 

//V 1 /^ //v^/v 2 ^ 2VTMVT + q(2 - VT) - q 2 (1 - VT)] 

<*o"o )s - (N Q N o)s = [2 _ Q(1 _Vf)]4 ' (B22) 



(ns*2*$)s = - " [ ^r^: v ^;; v± ~ v ^ (B23) 



2T 2 a[2VT + 2a(l - VT) - a 2 (l - VT)] 
[2 - a(l - VT)]S 

(TV 4T /V 2T ]V 21 )s = 0. (B24) 
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For energies fi < E < eV we have that (N^) B , (N*) B , {N$N#)b, (N^N^ l ) B , (N^N^N^) B and (n£n$N* 1 )b 
are equal to the case with two additional reservoirs. The other expectation values are 



T 
2' 
= T, 

2T 2 (l-a)sin 2 (0/2) 



<W>»- [2 _ a(1 _^)] 



(NkNg)B = (Nq Nq) b = 



2T 2 (1 - a) 



[2 - a(l - v^)] 2 



(B25) 
(B26) 
(B27) 

(B28) 



APPENDIX C: PROBABILITY DISTRIBUTIONS 



In order to calculate the various probabilities needed to evaluate the CH inequality. Eq. 1)12(1 and ((20(1 . it is necessary 
to solve the integrals of Eq. ((14(1 , where the different characteristic functions are given in Eq. JED) and Eq. |B2)l . As 
we mentioned in Section 1 1 1 1 1 explicit expressions of the probability distributions in terms of sums are lengthy and 
complicated for practical calculations. It is possible, nevertheless, to express the result for the various probabilities 
in a quite simple fashion, which makes them manageable for computational analysis. The point is to realize that the 
characteristic functions are nothing but polynomial functions on the variables e lAi of different degrees. The effect of 
each integral of Eq. 1(14(1. together with its accompanying complex exponential e lXi< ^ i /(2n), is simply to select the 
coefficient of the characteristic function polynomial which corresponds to the power equal to Qi. Consequently, the 
result of the integrals for the various probabilities can be expressed in terms of Qi order derivatives of the characteristic 
function, as we show below. Note that, whenever both analyzers are present, we will choose without loss of generality: 
a =1 and a' =\. For single terminal probability distributions we will also set a =f for lead 1 and a' =f for lead 2, 
although resulting expressions will not depend either on the direction of the spin or on the angle of the analyzer. 

For the single terminal probability distribution with analyzer we have 



i 

27 



dA 1T e !Al ^T x ( AlT ) 



r x(Ai T ) 



Q 



Qi 



(CI) 



where x(Ai|) can be extracted from Eq. ((Bl(l making A 2 j = and using Eq. 1(18(1 . Since fx = ^f, we have that 
= 4f , being M the total number of emitted particles per lead or per spin. The expectation values needed in 
Eq. ((Bl(l above and below energy /i are given in Appendix [B] The single terminal probability distribution in the 
absence of analyzer is 



P(Qi 



1 

2^ 



dA ie jAiyi x (Ai) 



1 



L A(Ai 



iA i->0 



(C2) 



where the characteristic function can be extracted now from Eq. (|B2|) . setting again A 2 f = 0. We can get similarly 
the expressions for P 2 (<32t) an d P{Q2)- 

The joint probability distribution when both analyzers are present gives 



which only depends on the angle O 
have 



1 



(2^) 2 
1 



dAite^Qn / rfA 2T e lA2 T^T x ( AlT ,A 2T ) 

J — 7T 



Qi T !Q 2T ! d(e 4A iT)QiTd(e lA 2T)Q2 T 



e ' A l T e iA 2 T ^ 



(C3) 



?i ± 62) /2, as we showed in Section ITTT1 When there is only one analyzer we 



and 



1 



IT: 1 



p-' e2 (QuQ2i) 



(2^ 



(2n)i 



dAi T e lAlTQlT 



rfAie 



iXiC 



dA 2 e 4A2< 3 2 X (Ai T ,A 2 ) 



dA 2T e lA2 T%T x ( Al ,A 2T ) = 



! x(Ait, A 2 ) 



Qi T !Q 2 ! d(e iX ^)Q^d(e iX i) c 



2T x(Ai, A 2T ) 



(C4) 



e * A lT , e iA 2^0 



Qi!Q 2T ! d(e lA i)Qid(e a2 T)Q2 T 



.(C5) 
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However, this two last expressions are not strictly needed since one can use the relations in Eq. (|19f) . Again, all the 
expectation values which are needed for evaluating these pr obabilities are given in Appendix iBl 
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